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Phase Diagram of the BCC S — \ Heisenberg Antiferromagnet 
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We use linked-cluster series expansions, both at T = and high temperature, to analyse the 
phase structure of the spin-i Heisenberg antiferromagnet with competing first and second- neighbor 
interactions on the 3-dimensional body-centred-cubic lattice. At zero temperature we find a first- 
order quantum phase transition at J2/J1 — 0.705 ± 0.005 between AFi (Neel) and AF2 ordered 



phases. The high temperature series yield quite accurate estimates of the bounding critical line 
for the AFi phase, and an apparent critical line for the AF2 phase, with a bicritical point at 



J1/J2 ~ 0.71, kT/Ji — 0.34. The possibility that this latter transition is first-order cannot be 
excluded. 

I. INTRODUCTION 



> 

The occurrence of competing exchange interactions in magnetic materials can give rise to a rich variety of magnetic 
ordered states, and of phase transitions between them. Studies of such phenomena, within the "molecular" or "mean- 
field" approximation go back half a centuryi or more. It is perhaps surprising that open questions remain, but, at 
£C) ' least for quantum models, this is the case. 



We study in this paper, the spin-i Heisenberg antiferromagnet on the body-centered-cubic (bcc) lattice, with first- 



' and second-neighbor interactions. The Hamiltonian is 

s 

^: ff = j 1 ^ (1) s i -s J +j 2 ^ (2) s,-s j a) 

£ 1 (ij) («) 

o' 

where the summations are over first and second- neighbor pairs respectively and the exchange constants J\,Ji > 0. 

. For our purposes we need to divide the structure into four interpenetrating sublattices, which we denote A\, A2, B\, 
• 1— 1 . 

B2 and illustrate in Fig. 1(a). Each of these sublattices has a face-centered-cubic (fee) structure. The J\ interactions 
couple the A and B sublattices, while the interactions J2 couple A\ to A2, and B\ to B2 only. 

The first question to ask concerns the ground state, at T — 0. For classical vector spins it is easy to show that this 
will be the Neel or AFi phase for J2/J1 < 2/3, and the AF2 phase for J2/J1 > 2/3. In the AFi phase all A spin point 
in the direction of an arbitrary unit vector n while B spins point in the opposite direction — n. In the AF2 phase each 
of the A and B sublattice is itself Neel ordered. Thus we may choose A\ spins in direction n, A2 spins in direction 
— n, B\ spins in direction n' and B2 spins in direction — n', with n and n' completely independent. In the quantum 
case these will not be exact eigenstates but, as usual in quantum antiferromagnetism, will be modified by quantum 
fluctuations. Furthermore in the AF2 phase the direction n' will be locked to either ±n, leading to a discrete 2-fold 
Ising symmetry added to the 0(3) spin symmetry. 

This model has similarities to the so-called "Ji — J2 model" on the square lattice, which has been much studied 
in recent year o 2 ^ 4 . In that case there is now strong evidence for an intermediate "spin-liquid" phase between the 
Neel and AF2 (the "collinear phase"), with the possibility of an even richer structure. Motivated by this similarity, 
Schmidt et al£ have studied the bcc model at T = 0, using exact diagonalizations and linear spin- wave theory. They 
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FIG. 1: (a) The four sublattices of the bcc structure, (b) the AFi amd AF2 ordered phases. 



find a direct first-order quantum phase transition between AFi and AF 2 phase at J2I J\ — 0.7, with no intermediate 
phase. This can be understood in terms of the diminishing effect of quantum fluctuations in the 3-dimensional system. 

In Section 2 we use linked-cluster series expansions at T = to study ground state properties of the model. We 
obtain rather precise results for both the ground state energy and order parameters, confirming the scenario of Ref. 
P| and obtaining a somewhat more precise estimate of the quantum phase transition point. 

A second interesting aspect of the model is the phase structure at finite temperatures. There will be transition lines 
in the (T, J2/J1) plane where the low-temperature ordered phases meet the high-temperature disordered phase, and 
the locations and nature of these transition lines are not well known. The earliest attempt to address these questions, 
beyond mean-field theory, was the 6th-order high-temperature series work of Pirnie et al&. They obtained series 
for the appropriate staggered susceptibilities for both AFi and AF2 ordering for the J\ — J2 bcc lattice (for general 
spin), and estimated critical temperatures via Dlog Pade approximants, in the usual way. They were unable to obtain 
precise values for the critical exponent 7. More recently Pan£ has computed an 8th-order staggered susceptibility 
series for the Neel phase for the case J2 = 0, and obtained the critical temperature as kT c j J\ = 1.384 ± 0.005. (This 
is actually a thermodynamic perturbation expansion, which can generate a high-temperature series but can also used 
to study properties in the ordered phase). 

In Section 3 we report on the derivation and analysis of new extended high temperature series for this system. In 
particular we have computed staggered susceptibilities for the full J\ — J2 model to 10th order in the AFi phase and 
9th order in the AF2 phase, adding 4 and 3 terms respectively to the work of Ref. jfj. These longer series allow 
reasonably precise estimates to be made of both critical temperature and exponents. There is one caveat. There 
are arguments that the AF2 to paramagnetic transition is first-order. If that is indeed the case then our series are 
presumably seeing a spinodal line within the ordered phase, and the actual transition temperature is higher. 

The arguments for a first-order transition come from two sources. Firstly for the Ising version of this system both 
Monte Carlo simulations^ and a combination of high and low-temperature series 9 support, fairly conclusively, the 
existence of a "fluctuation-induced" first-order transition. Secondly, a renormalization group treatments*^ of the 
appropriate n — 6 component Landau- Ginzburg- Wilson Hamiltonian for the AF2 transition has no stable fixed point, 
suggesting a first-order transition. However this analysis is based on an e-expansion about d = 4 and, in our view, 
while persuasive, is not conclusive. 
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FIG. 2: Ground state energy versus J2/J1 from Ising expansions for the AFi and AF2 phases. The crossing point at 0.705±0.005 
identifies the first-order quantum phase transition. Uncertainties in the series extrapolation are no larger than the symbols. 
Lines joining symbols are guides to the eye. 

II. EXPANSIONS AT T = 

Our approach is to use high-order linked-cluster expansions^ where the Hamiltonian is written as 

H = H + xV (2) 

where Ho is the Ising Hamiltonian 

H = J, J2 SfS] + J 2 J2 ^ SfSI (3) 

(ij) (ij) 

and 

v = \ J i E (1) ( s ? s 7 + S ^ S T) + \ J * E <2) ( s ? s i + s ^ s i) ( 4 ) 

(ij) (ij) 

is treated as a perturbation. Quantities are expanded as power series in x, and then evaluated at the isotropic limit 
x = 1, by Pade and differential approximants. We have used this technique successfully in many previous studies of 
quantum antiferromagnets^Siii. In the present work we are able to compute series through x a in both Neel and AF2 
region. 

Figure 2 shows the ground state energy versus exchange ratio J2/J1, the two branches corresponding to the two 
phases which can occur. The results „ ve ry to those of Ref. Q obta.oeo fr o ra extrapolate of exact 

diagonalization results. The two branches clearly cross, confirming that the ground state transition between AFi and 
AF2 phases is first-order. Our estimated transition point is J2/J1 — 0.705 ± 0.005. 
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FIG. 3: Staggered magnetizations (order-parameters) and correlator discriminants (Eq. |5J versus J2/J1, in both the AFi and 
AF2 phases. ACi , AC2 refer to first and second neighbour correlators respectively. In the AF2 phase the superscripts ++ and 
+ — refer to first neighbour correlators for sites with like spins and unlike spins. The thin vertical line denotes the transition. 
The lines joining symbols are guides to the eye. 



We have also computed series for two other physical quantities: staggered magnetization in each phase M, and the 
correlator discriminant 

&C=\3(S;S5)-(3 i .S j )\ (5) 

which is a measure of the breaking of spin rotation symmetry. These results are shown in Figure 3. We look first 
at the magnetization. In the Neel phase this starts from ~ 0.435 and decreases monotonically to approx. 0.34 at 
the transition point. In the AF2 phase the curve start from 0.424, the result for the simple cubic lattice^, and 
remains quite flat, with a downward curvature near the transition point. Our results show a small upward curve near 
J2/J1 ~ 1, but this is most likely a artifact. Our results are similar to, but considerably more precise than, those of 
ref. fl. 

Results for the correlators are new. Rather than show all correlators, we show in Figure 3 only the quantities AC in 
Eq. [5] for first and second- neighbor pairs. A nonzero value indicates spontaneous breaking of spin rotation symmetry. 
Of course our expansion starts from a symmetry-broken state, and is thus biased. However an identical calculation 
for the Ji — J2 square lattice^ clearly shows AC — > at the boundary between the Neel and disordered phase. The 
fact that this does not happen here is clear and strong evidence that the Neel phase does not vanish continuously, 
but undergoes a first-order transition to the AF 2 phase. The second-neighbor discriminant AC2 in the AF 2 phase 
is particularly interesting, showing a small but clearly discernable increase on approaching the transition. This is 
unexpected. 

In summary our T = series confirm the occurrence of two phase, AFi and AF 2 with a first-order quantum phase 
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FIG. 4: Phase diagram of the BCC Ji — J2 antiferromagnet. The solid line is the AFi-paramagnetic critical line, with universal 
n = 3 exponents. The dashed line is the AF2-paramagnetic transition line, as determined from our series. The short vertical 
dotted line is the first-order AF1-AF2 transition line, extending from the bicritical point (or crtical end-point) to the quantum 
phase transition point on the T — axis. 

transition at Jij J\ = 0.705 ± 0.005. We now turn to the finite temperature phase structure. 



The standard approach to determining finite temperature critical points and exponents for quantum (and classical) 
models is via high-temperature expansions. For the bcc(l,2) antiferromagnet this was first done in Ref. 6]. One 
includes an appropriate staggered field F in the Hamiltonian 



where rji = ±1 on the different sublattices, reflecting the type of order expected at low temperatures, expands the 
partition function 



III. HIGH TEMPERTURE EXPANSIONS 




(6) 




OO 



(7) 



n=0 



and compute the appropriate staggered suscetibility 



x = - 




dF 2 



(8) 



in the form 



DC 




(9) 
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where K — 0J\, a — J2/J1 and a constant has been factored out. Ref. Q gives the a n (a) for both AFi and AF2 
phases through n = 6. We have extended these series by 4 terms in the AFi phase and 3 terms in the AF2 phase. The 
coefficient are given in Table I. Our results agree fully with previous shorter series&£ for the AFi phase, but there is 
a discrepancy at 6th order for the AF2 series. We believe that in the second last row of the last table in Appendix 
II of Ref. |6J the entry 994944 should be replaced by -2230656. With this substitution we find agreement for both 
S = h and S — 1 (for which we have also derived short series). 

The series have been analysed via standard Pade approximant methods 14 . Table II shows the critical temperature 
and exponent estimates obtained from a direct analysis of -^\nx(K). As is apparent, for the smaller value of 
a = J2/J1, there is a consistent pole and an exponent around 1.4, which is the value expected for the Heisenberg 
(n = 3) universality class. As a increases the series become less regular. This is due to a pole on the negative 
real axis lying closer to the origin ("ferromagnetic singularity"). Nevertheless the direct analysis is consistent with 
7 ~ 1.4 along the entire line. Assuming this, it is possible to increase the precision in K c . The resulting critical 
line is shown in Fig. 3. In the AF2 phase the series are less regular and consequently the analysis is less precise. In 
Table III we show some of the raw analysis results. As is evident the higher-order Pades show a consistent physical 
singularity down to at least J2/J1 = 0.8, although there is again an interfering singularity on the negative real axis. 
The exponent estimates show a clear decreasing trend from J\/ J2 — (while corresponds to the Neel phase of the 
simple cubic lattice, for which we expect a second-order transition with the universal exponent 1.40) for increasing 
J1/J2. For example, for J1/J2 = 1.25 (J2/J1 = 0.8), the estimated 7 < 1. Now, if the AF2-paramagnetic phase 
boundary is second-order, we would expect critical exponents appropriate to an n = 6 component order parameter. 
Thus the series will be affected by a crossover from n = 3 to n = 6, presumably followed by a second crossover to the 
bicritical point. In a relatively short series this can give rise to varying critical exponents, as observed here. On the 
other hand there are arguments, mentioned previously, for this transition to be first-order. If that is the case then 
the transition will lie above our curve in Fig. 4, and the apparent divergence locates a spinodal line within the order 
phase. The "bicritical point" is then, in fact, a "critical-end point"—. 



IV. CONCLUSIONS 

We have used a combination of series expansions at T = and conventional high-temperature series to elucidate 
the phase diagram of the quantum S — | Heisenberg antiferromagnet on the BCC lattice with nearest and next- 
nearest-neighbor interactions. The ground state properties of this model have not, to our knowledge, been previously 
studied by series methods. We identify a first-order quantum phase transition at J2/J1 — 0.705 ± 0.005 between 
the two possible types of antiferromagnetic order. We have extended previous high temperature series by 4 and 
3 terms, respectively in the AFi and AF2 phases. These longer series allow a fairly precise estimate of the AFi- 
paramagnetic critical line, and confirm n = 3 universality along this line. The nature of the AF2-paramagnetic line 
remains enigmatic. We find a consistent line of poles, which may represent a true critical line. However we cannot 
exclude the possibility of a first-order transition at somewhat higher temperatures. 

In the corresponding Ising case the AF2 -paramagnetic transition was shown to be first-order by a careful analysis 
using both high and low temperature series^. The former approach cannot be used here since low temperature series 
cannot be obtained. It would be of some interest to study this system using quantum Monte Carlo methods, although 
a "sign problem" would be expected. A thorough Monte Carlo study for the classical Heisenberg system would also 
appear very worthwhile. To our knowledge no such studies have been attempted. 
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TABLE I: Susceptibility polynomials a n (a) (Eq. for the spin-i bcc(l,2) Heisenberg antiferromagnet. To avoid fractions, the 
values of 2 n+1 (n + l)\a n (a) are given. 

AFi (Neel) phase 

n=l 16, -12 

n=2 160, -288, 72 

n=3 2048, -5424, 4032, -528 

n=4 31096, -110640, 132168, -55680, 4950 

n=5 557456, -2470992, 4070976, -2947488, 825120, -56196 

n=6 11495460, -61270752, 127951356, -131155632, 65324256, -13458816, 693273 

n=7 269007424, -1670678976, 4241825472, -5589739344, 4034541024, -1489565472, 238760256, -9630816 
n=8 7026835032, -49911466560, 148837552176, -240765473760, 226149075384, -123139222944, 35480202336, 
-4568142528, 156881934 

n=9 202835216096, -1620119828784, 5556372919872, -10631250963456, 12358650036192, -8857806164208, 
3806668562496, -888255818640, 94585942080, -2810097960 
n=10 6406925312668, -56883245890656, 220319075124360, -487248245389392, 674224356164040, -602298324975312, 
344178723975168, -120654960988368, 23412924414240, -2116355382240, 52557775149 

AF2 phase 

n=l 0, 12 

n=2 -40, 0, 84 

n=3 -64, -240, 0, 720 

n=4 1320, -528, -5112, 0, 7422 

n=5 9968, 3408, -11472, -50784, 0, 92412 

n=6 -131916, 13200, 197604, -133248, -1262928, 0, 1323657 

n=7 -2026752, 6866304, 2195328, 12122160, -3610560, -18223392, 0, 21144864 

n=8 556840, 68279328, -335099184, 70461024, -2954040, -61043232, -432565248, 0, 375695526 

n=9 359570256, 205680, -3444654528, 10167197376, 740209344, 8510348784, -1552929696, -7905850224, 0, 7456498848 



TABLE II: Estimates of the critical temperature K c and exponent 7 (in brackets) from Dlog Pade approximants to the Neel 
phase staggered susceptibility series. 



J2/J1 0.0 0.2 0.4 0.6 

[3/4] 0.7233(1.391) 0.8881(1.427) 1.1614(1.418) 1.6710(1.117) 

[4/4] 0.7248(1.410) 0.8839(1.391) 1.1634(1.429) 1.7446(1.314) 

[3/5] 0.7251(1.414) 0.8844(1.396) 1.1635(1.430) 1.7706(1.416) 

[5/4] 0.7281(1.459) 0.8896(1.454) 1.1704(1.478) 

[4/5] 0.7199(1.371) 0.8859(1.410) 1.1606(1.415) 1.6534(1.087) 

negative pole -1.0 -0.8 -0.68 -0.6 
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TABLE III: Estimates of the critical temperature K c = J^jkT^ and exponent 7 (in brackets) from Dlog Pade approximants to 
the AF2 phase staggered susceptibility series. 



■hi 32 


0.0 


0.5 


1.0 


1.25 


[3/3] 


1.0218(1.232) 


1.0698(1.299) 






[2/4] 


1.0318(1.287) 


1.0709(1.304) 






[3/4] 


1.0540(1.404) 


1.0847(1.367) 


1.221(1.232) 


1.387(0.966) 


[4/4] 


1.0627(1.461) 


1.0827(1.357) 


1.233(1.282) 


1.369(0.920) 


[3/5] 


1.0639(1.472) 


1.0828(1.358) 


1.234(1.288) 


1.370(0.924) 


negative pole 


-1.5 


-0.7 


-0.5 


-0.4 



